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PACS numbers: 

A major motivation for cold atom experiments 
is the search for quantum ground states such as 
antiferromagnets and d-wave superfluids. The 
primary obstacle to this task is the difficulty of 
cooling to sufficiently low temperatures. We pro- 
pose a way to achieve very low temperatures 
and entropies O.OSfcs per particle) by trapping 
fermions in a corral formed from another species 
of atoms. The Fermi system can then be used 
as a heat sink, or it can be adiabatically evolved 
into other desired states. In particular, we sug- 
gest methods for generating antiferromagnetism 
using this technique. 

The last decade has seen astonishing progress in ul- 
tracold atom technology. In particular, atomic gases in 
optical lattices can now be imaged and manipulated with 
amazing controllability and fidelity [1^4], in a way that 
parallels the development of atomic-resolution scanning 
tunneling microscopy in the 1980s. 

Recent advances in the Harvard group [5] have suc- 
ceeded in creating a large bosonic Mott insulator (MI) 
with one atom per site. The core of this MI can be said 
to have zero entropy density. Naively, one might think 
that the MI core has zero temperature. However, this 
is not true, because this is not the ground state of the 
Hamiltonian with the particular trap potential. Thus, if 
the lattice is lowered adiabatically in an attempt to ac- 
cess interesting quantum ground states, the populations 
in the center of the trap will readjust, causing heating. 
A more fundamental limitation is that entropy from the 
"fuzzy" cloud boundary, where the occupation number 
is randomly or 1, will leak inwards and contaminate 
the MI. To minimize this effect it is desirable to engi- 
neer a system in which the boundary region is as small 
as possible. Recent advances in laser beam shaping using 
digital micromirror devices have been able to create flat- 
top beam profiles with less than 1% RMS error [6], which 
may be used to create "square wells" for atomic gases. 
However, optical beam shaping is inherently diffraction- 
limited and would inevitably produce a boundary region 
a few sites thick, where the number is uncertain and there 
is unwanted entropy. 

Preparation of zero-entropy states: Here we pro- 
pose a method that can, in principle, create pure states 
that are also eigenstates of the Hamiltonian. The key 
idea is to form a corral using atoms instead of light. This 
allows the corral potential to be controlled with single- 




FIG. 1: Fermions in a bosonic corral: Large green 
spheres represent bosons and small blue spheres represent 
fermions. The grid represents the optical lattice. The 
fermions are trapped in the corral and are unable to tunnel 
out. In the absence of vacancies (such as the one purposely in- 
cluded in this illustration) the fermions form a zero-entropy, 
zero-temperature ground state of an isolated system. This 
is like using a cookie cutter to isolate a clean, homogeneous 
piece of a slab of dough. 

site resolution, so that the boundary region is essentially 
eliminated. 

This idea may be implemented in ID, 2D, or 3D; we 
will focus on the 2D case. A single species of fermions 
is loaded into a trap, and an optical lattice is gradually 
ramped up [27]. The resulting 2D cloud consists of a band 
insulator (BI) core surrounded by a Fermi liquid (FL) 
shell. A bosonic Mott insulator is loaded into the same 
spatial lattice and patterned into a hollow 2D shape [2-4] . 
Experimental parameters are chosen such that the bosons 
experience a deep lattice {V^'^^^ > 20i?^), so that boson 
tunneling is negligible and the bosons stay in place for 
the remainder of the experiment [28] . The boson-fermion 
interaction L^bf is now ramped up to a large repulsive 
value. The bosons form a potential barrier of height C/bf 
that is impermeable to the fermions. The N fermions 
within this "corral" now constitute an isolated quantum 
system with exactly one fermion per site (see Fig. 1). 
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Because the state is uniquely specified, it has zero en- 
tropy. Furthermore, this state is the ground state of 
the Hamiltonian in the subspace of total particle num- 
ber TV. This is consistent with a Boltzmann distribution 
with temperature T = 0, and we do indeed have a zero- 
temperature system in its equilibrium ground state. 

Corral shape and preparation: The corral of 
bosons may be prepared using the method of Ref. [4], 
where a tightly focused laser beam, together with a mi- 
crowave field, was used to flip the hyperfine spin of indi- 
vidual atoms in a singly-occupied bosonic Mott insulator 
with sub-diffraction-limited resolution well below the lat- 
tice spacing, and a push-out laser pulse was then applied 
to remove the flipped (or unflipped) atoms. 

Alternatively, as in Refs. [2, 3], one may begin with 
a large bosonic system consisting of a doubly-occupied 
Mott insulating core surrounded by a singly-occupied 
Mott insulating ring, and then illuminate the atoms to 
eject pairs of atoms via light-assisted collisions. This re- 
sults in a core of empty sites surrounded by a MI ring, 
which can then be used as the corral. Although this tech- 
nique provides less control over the exact shape of the 
corral, it has the advantage that it may be used to pro- 
duce a spherical corral which can then be used to make 
a 3D antiferromagnet (AF), which is much more robust 
than its 2D analog. 

An appealing feature of our proposal is that the corral 
does not have to be in the shape of a perfect square, cir- 
cle, or sphere. The only requirement is that the bosons 
form a multiply connected region with a continuous wall 
enclosing some empty sites. The digital (in this case, bi- 
nary) nature of the boson occupation number leads us 
to a zero-entropy fermionic state, regardless of irregular- 
ities, as illustrated in Fig. 1. 

In order to confine the Fermi band insulator, the cor- 
ral potential height f/bf should be chosen to be greater 
than the bandwidth. Also, to prevent fermion loss by 
quantum tunneling through the walls, and to minimize 
the probability of gaps and leakage, the corral should be 
made several bosons thick. 

Discarding excess fermions: Although the interior 
of the corral is already an isolated system, it is neverthe- 
less a good idea to get rid of the excess fermions on and 
outside the corral to prevent them from interfering with 
later stages of the experiment. Figure 2a shows a slice of 
Fig. 1 in the x-direction. Up to this point the fermions 
are still under harmonic confinement. Now, the fermionic 
trap potential in the {x,y) plane is inverted (while main- 
taining z-confincmcnt). That is, it is changed to 
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Then, the exterior fermions accelerate outward and 
rapidly exit the system, whereas the interior fermions 
remain confined by a potential barrier that is too thick 
for them to tunnel through. (It may be advantageous to 



6^ 



FIG. 2: Discarding excess fermions by trap inversion: 

(a) First, the Fermi insulator is formed in a harmonic trap, the 
bosonic corral is created, and the boson-fermion interaction 
is ramped up. Black curves show the combination of the trap 
potential, y^^^^ ^ and the Hartree potential, f/bf?i6, seen by the 
fermions. (b) The trap potential is now "turned upside down" 
to "pour away" unwanted fermions. (c) Finally, the inverted 
trap potential is turned off to give a flat-bottomed infinite 
square well. The lattice and ^-confinement are maintained 
throughout this process. The process is akin to throwing away 
unused dough after cutting out a cookie of the desired shape. 



lower the optical lattice depth to facilitate this separa- 
tion.) After a suitable time, the trapping frequencies are 
reset to zero in the x and y directions, i.e.. 



trap / 
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This leads to the ideal situation of fermions in a flat- 
bottomed square well with infinitely high walls. We 
have performed simulations of this trap-emptying proto- 
col (see supplementary information), which suggest that 
it is indeed realistic and feasible [29]. 

Adiabatic evolution to interesting ground 
states: An = T = system can in principle be used as 
a heat sink to cool another system. However, the single- 
species Fermi insulator described above would be an in- 
efficient heat sink (due to its excitation gap, low heat 
capacity, low thermal conductivity, limited phase space); 
BECs and other fluids generally function much better as 
heat sinks. It is therefore most desirable to look for ways 
to adiabatically evolve the S — T — ^ system to an in- 
teresting ground state. In the remainder of this paper we 
present a way to attain an antiferromagnetic state. 

Protocol for obtaining a 2D antiferromagnet: 
Two species of fermions (f and J,) are loaded into a 
bosonic corral as previously described, such that every 
site of the square lattice is doubly occupied (1 \). (See 
supplement for a discussion.) The fermions experience an 
on-site, inter-species repulsion U ~ 6t (which may have 
been introduced at any stage in the protocol). The op- 
tical lattice is now gradually morphed into a lattice with 
twice the number of sites [30] as illustrated in Fig. 4. 
This may be done by ramping up a sublattice or by 
ramping down a superlattice. The feasibility of studying 
cold atom systems in double-well potentials is well es- 
tablished, having been demonstrated experimentally by 
several groups [11-15]. (During the preparation of this 
manuscript, we noted that such an idea has been pro- 
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FIG. 3: Splitting a singlet using a double-well poten- 
tial: Beginning witii a pair of fermions in a well of a deep opti- 
cal lattice, which is then adiabatically morphed into a double- 
well shape, one obtains a two-site singlet ^(It)^!^)^ ^ 

li>Jt)«). 




FIG. 4: Site-doubling route to antiferromagnetism: A 

square lattice is loaded with two fermions per site, forming 
a band insulator (a). Each site of the lattice is then split 
into a double well. The ground state consists of singlets in 
each double well (b). As the hopping between double wells 
(t') increases, antiferromagnetic correlations grow, and at a 
certain value of t' /t there is a quantum phase transition to 
AF order. 

posed very recently [16].) 

Starting from a fully occupied 2D Hubbard model at 
filling (n) = 2, the site-doubling procedure leads to a 2D 
Hubbard model at half-filling, (n) — 1. If this procedure 
is performed adiabatically, it will lead to the ground state 
of the final model, which is a singlet state with strong 
antiferromagnetic correlations. [31] 

The AF correlations will manifest themselves in the 
noise spectrum. For example, snapshots of either spin 
species obtained using an atomic gas microscope will 
show a checkerboard-like pattern. Alternatively, one may 
perform Bragg scattering to search for a peak at the an- 
tiferromagnetic wavevector, which is (tt, tt) in final recip- 
rocal lattice units. 

Comparison with previous suggestions for 
studying quantum ground states: Many cooling 
methods have been proposed and/or implemented so far. 
Some of these methods involve open systems with driv- 
ing and dissipation (laser cooling, Zeeman slowing, evap- 
orative cooling, and sympathetic cooling) [18, 19], while 
others operate on isolated systems (adiabatic expansion, 
"Pomeranchuk" cooling, and other adiabatic processes 
[20-22]). All of these methods involve starting with a 
system at a finite temperature Ti and entropy 5*1 , which 



is then evolved to a state of lower T and S. 

In our proposal, by "cordoning" off a portion of the 
system, we ensure that our initial state is already at T = 
and 5* = (in principle). Hence, we only need to worry 
about heating due to non-adiabaticity. If adiabaticity 
can be maintained, then we will arrive at a final state 
with T = and S* = 0. 

Estimates of deviations from ideal behavior: We 
finally attempt to estimate the entropy and temperature 
of the final state due to non-idealities at various stages 
of the protocol. 

(a) Imperfections in band insulator: The first 
stage of the process involves preparing a Fermi band in- 
sulator in the center of a trap (by evaporative or sympa- 
thetic cooling). Here, we benefit naturally from entropy 
distribution in an inhomogeneous system: the outer shell 
is a gapless Fermi liquid that can "soak up" entropy, 
so that the core, which is a gapped Fermi band insu- 
lator, has a lower entropy density than average. As 
an example of typical parameters, we consider ^^Rb in 
a 3D cubic lattice of spacing d — 680 nm with trap 
frequency ^Ixy = 50 Hz and lattice depth Vq — SEj^. 
In this case the nearest-neighbor hopping amplitude is 
thop/ks = 1.9 nK and the trap constant is a « 0.003ihop, 
such that V{r) = air/dY'. For a non-interacting Fermi 
gas with N — \{f particles and overall entropy per par- 
ticle S/N — 0.65fcB, the density and entropy profiles are 
shown in Fig. 5. Using a spherical corral of radius 25(i we 
can isolate a core subsystem with vacancy concentration 
0.008 and average entropy per particle 0.03/cb. While 
this entropy is not strictly zero, it does represent a 25- 
fold reduction from the starting average entropy, and it 
is certainly well below the critical entropy for antiferro- 
magnetic ordering in 3D [21, 22]. 
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FIG. 5: Density profile: The blue curve is the density pro- 
file of non-interacting fermions in a cubic lattice in a harmonic 
trap, according to the local density approximation. The sys- 
tem consists of a n = 1 band insulator surrounded by a Fermi 
liquid outer shell with < n < 1 that contains most of the 
entropy. The gray region represents a bosonic corral of inner 
radius 25ci and thickness 5d. 

For a 2D square lattice the entropy redistribution ef- 
fect is less dramatic (because the size of the Fermi liquid 
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region is amplified by the Jacobian, which is now only 
27ri? instead of AttR^). On the other hand, it should be 
easier to reach lower initial temperatures for a 2D system 
by evaporative cooling, because one can evaporate away 
all the atoms that are not in the plane of interest. 



FIG. 6: Heisenberg model on dimer lattice: Thick 
bonds represent exchange coupling J and thin bonds repre- 
sent coupling J' < J. 

(b) Non-adiabaticity in site-doubling process: 

We now study the possibility of adiabaticity/non- 
adiabaticity in the process shown in Fig. 4. As we go 
from isolated dimers (Fig. 4b) towards a rectangular lat- 
tice (Fig. 4c), the system remains half-filled with one 
fermion per site on average, and the "charge gap" re- 
mains robust. Thus, for the purpose of modeling the 
process we focus on the spin degrees of freedom. We 
consider a SU(2) = 1/2 Heisenberg antiferromagnet on 
a lattice of dimers, with coupling J within each dimer 
and J' between dimers, 

^^Hcis = + ^ JijSi ■ Sj, (3) 

(ij) 

where (ij) are pairs of nearest-neighbor sites and J.y- are 
couplings as described above. [32] We start with J — 1 
and J' = (isolated dimers) . The ground state is a spin- 
gapped singlet state, with singlet "valence bonds" living 
on strong bonds J^, which we shall call a "valence bond 
solid" (VBS) for want of a better term, 

m^o))^ (g) ;^(it)j;),-i;)jt),). 

ijeCi 

We represent it in the Sz eigenbasis, as a vector with 
2''^='*°= coefficients. Now the hoppings = Jy are 
ramped from to 1 as a function of time t over a to- 
tal time iramp- This goes from the dimer limit to the 
isotropic square lattice limit. We evolve the wavefunc- 
tion by solving the time-dependent Schrodinger equation 
^ |i/)(i)) = -iH{t) iV'(i)) with the Verlet method, where 
H and \^) are time-dependent (see supplementary infor- 
mation). As the lattice becomes more isotropic, anti- 
ferromagnetic correlations grow, which is an indication 
of the quantum phase transition to an antiferromagnetic 
state in a bulk system. 

The staggered structure factor Sq — (Afq^), where 
AIq = ^-(expzQ • 'R.i)Sf and Q = (7r,7r), indicates the 



size of antiferromagnetic fluctuations. With this defini- 
tion the "VBS" state has Sq = 7Vsites/2, whereas a clas- 
sical Neel antiferromagnet with a Neel vector Mq in the 
z direction would have Sq = iVsites^/4. Averaging over 
all directions of Mq on the surface of a sphere reduces 
this to Sq = 7Vsites^/12. Quantum fluctuations further 
affect this number. 

We also calculate the mean energy Savg — (V'l^l'/') 
and the energy variance i?wid^ — ("01 H'^ IV") ~ (V'l l'/')^- 
i?wid is a measure of the admixture of excited states pro- 
duced by non-adiabaticity; roughly speaking, it is also 
an indication of the eventual temperature due to heating 
due to non-adiabaticity (although we do not model the 
thermalization process here). 

Figure 7 shows -Eavg: ^-wid, and Sq as a function of 
time t for various ramp rates 1/iramp- For sufficiently 
slow ramp times the behavior of these quantities shows 
that the evolution is adiabatic (see caption). This al- 
lows us to estimate the fldelity (the squared overlap of 
the wavefunction with the instantaneous ground state) as 

f {tramp, J') = |(V'(30, J') |V'(^ramp, J')> 1^7 showu in the 

last panel of Fig. 7, which gives further evidence for our 
results. 

For a hopping amplitude thop = 2 nK and a Hub- 
bard repulsion U = 8<hop, the exchange coupling scale 
is J ~ 4t^Qp/f7 — 1 nK, corresponding to a time scale of 
1/21 Hz = 50 ms. Thus, the ramp time should be at least 
500 ms. This is short enough to permit measurements in 
the antiferromagnetic state to be made before the state 
is destroyed by atom losses and/or technical heating. 

Adiabaticity is relatively easy to achieve in the disor- 
dered phase, where there is a spin gap. In a bulk system 
the spin gap decreases to zero at the QCP, and the ex- 
citation density scales as a power law of the ramp time, 
governed by appropriate critical exponents [24-26] . How- 
ever, in a finite system of linear size L, the gap at the 
QCP scales as (because the dynamic critical expo- 
nent is z = 1). Therefore the ramp time necessary for 
adiabaticity scales only as L. This is a likely explanation 
for the observations of Ref. [16]). For sufficiently large 
ramp times the excitation density is exponentially small 
according to Landau-Zener theory. This is an encourag- 
ing observation. 

Conclusion: In conclusion, we have proposed a very 
promising method for obtaining quantum ground states. 
The experimental challenge is to find a system where the 
lattices and interactions can be tuned appropriately. 

We gratefully acknowledge support from DARPA grant 
no. 60025344 under the Optical Lattice Emulator (OLE) 
program. We are grateful to William Cole and Eric 
Duchon for helpful discussions. 
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FIG. 7: Time evolution of various quantities: Mean energy i?avg, energy uncertainty iSwid, staggered structure factor 
Sq, and fidelity / on a 6 x 4 lattice as a function of interdimer coupling J' during the AF-generating protocol for various ramp 
times tramp- For iramp ^ ^/J (suddcu turn-on), the wavefunction has no chance to evolve, so the AF structure factor Sq stays 
constant and the energy uncertainty grows quickly. For tramp 10/J (adiabatic ramp-up), Sq grows indicating the growth 
of AF correlations across the quantum critical point, iSavg tracks the ground state energy, and i?wid remains small, indicating 
high fidelity (negligible admixture of excited states). The black dashed line in (c) indicates the structure factor for a classical 
Neel antiferromagnet Sq = N^itas/^'^- The unit of energy is the intra-dimer coupling J. For a comparison between diiferent 
lattice sizes, see supplement. 
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Removal of excess fermions: details 



Naively, one might think that momentary application 
of a steep inverse harmonic potential would quickly cause 
the exterior fermions to fly away in all directions. How- 
ever, when the fermions have accelerated to the Brillouin 
zone edge (with crystal momentum k = Tr/d) they Bragg- 
scatter to the opposite edge (— tt/o?), the group velocity 
changes sign, and the fermions start moving inwards in- 
stead (this is the phenomenon of Bloch oscillations) . 

A more suitable approach is to apply the anti-trapping 
potential for a short duration to impart an outward crys- 
tal momentum to the exterior fermions (noting that the 
maximum group velocity occurs in the middle of the dis- 
persion relation) , and then turn off the potential to allow 
the fermion cloud to expand ballistically. 

To confirm the feasibility of this idea it suffices to 
study the problem in ID. Figure 1 shows a simulation 
of the dynamics of 30 non-interacting spinless fermions 
in a 150-site tight-binding chain. (All energies are quoted 
in units of thop, and all lengths are in units of the lattice 
spacing d. In a realistic experiment ihop ~ (25 ms)~^ 
and d ~ 600 nm.) The system is initialized in the 
ground state of a harmonic trap potential V{r) = a{0)r^ 
where a(0) = -1-0.02. At time i = 0+ we turn on a 
corral potential consisting of two potential barriers of 
height 5 and width 4. The trap constant a{t) is var- 
ied according to the protocol shown in the inset. We 
see that the final state does indeed consist of a n = 1 
band insulator for 65 < a; < 82, with negligible den- 
sity immediately outside the corral. See ancillary file 
(TrapEmptyingAnimation.mov) for an animation of the 
evolution of the potential and density profile. 

In real experiments the optical lattice beam has a fi- 
nite waist. Therefore, once the atoms have travelled a 
sufficient distance from the trap center, they no longer 
experience Bloch oscillations and they should quickly be 
lost from the system. 

There are other possible ways of removing excess 
fermions, such as using appropriate radiofrequency or 
laser pulses to "blast away" fermions at corral sites. 




FIG. 1: Potential V{x) (top panel) and density profile p{x) 
(bottom) for 30 fermions in a 150-site tight-binding model 
with a time-dependent potential. Dashed curves show the 
potential and ground state density at time t = 0. Solid curves 
show the potential and density a,t t — 40, when the excess 
fermions have moved far away from the corral. 



Evolution into antiferromagnetic state: details 

We used the Verlet method for time evolution, in which 
the wavefunction is split into real and imaginary parts, 
xp = u + iv, and evolved according to the second-order 
difference scheme 



u(t + l5t) = u(t) 
v(t + 5t) = v(i) 
u{t + 5t) = u(i -\ 



-iH(t 
- iH(t 
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\5t) 

\5t)^r{t + 5t). (1) 



The timestep 5t was made small enough so that the os- 
cillations in the norm of the wavefunction, {ip \'ip), were 
smaller than 5%. This is shown in Fig. 2. 

Figure 3 compares the time evolution of the observ- 
ables i^avg, -E-wid, Sq, and / for various system sizes. We 
see that for a ramp time ^lamp — 10/ J the energy uncer- 
tainty i?wid is very close to zero and the fidelity / is very 
close to unity, indicating adiabaticity. 



2 



1.10r 

1.05 



im 1-00 




0.95 



0.03 0.1 0.3 1 3 10 30 
tJ (log scale) 

FIG. 2: Time evolution of the norm of the wavefunction on 
a 6x4 lattice for various ramp times tramp = 0.25/ J to 30/ J. 



ized two-component Fermi liquid, turning kinetic energy 
into heat in the process. Due to the above difficulties, 
we propose using a two-component band insulator as a 
starting state and performing a site-doubling procedure 
as described in the text. 

In the text we considered site doubling using a double- 
well geometry. A possible alternative involves morphing 
from a square lattice with a.dxd unit cell to a 45° x ^ 
lattice, as illustrated in Fig. 4. (Here we assume that the 
new lattice only applies to the fermions.) In this case, the 
number of sites in the new lattice is not exactly twice 
the original number of sites due to edge effects, so the 
final state is a hole-doped antiferromagnet. On the other 
hand, this geometry might be less susceptible to leakage 
(fermions slipping past bosons). 

Other methods for analyzing non-adiabaticity 



Other methods for adiabatic evolution to an AF 

In order to obtain interesting states of fermions one 
has to have at least two species of fermions. The sim- 
plest version of our proposal produces only one species 
of fermions. The easiest way to generate two species of 
fermions is to apply a 7r/2 pulse to evolve the fermions 
from the original hyperfine state |1) to into a superpo- 
sition ;^(|1) + 1 2)), and then ramp down the fermionic 
lattice depth (V/) to allow tunneling between sites. How- 
ever, this approach involves a Hamiltonian with a mas- 
sively degenerate ground state (now that we are in the 
Fock subspace with A^i = iV2 = NsUcs/'^), from which 
adiabatic evolution is useless. If instead we begin with 
fermions in a shallow lattice permitting intersite tunnel- 
ing, then the coherent superposition ^(|1) + 1 2)) is not 
the ground state of the Hamiltonian, but rather, it is a 
fully polarized state (polarized along the a;-axis in pseu- 
dospin space) , which will eventually relax to an unpolar- 



Apart from analytic theories (Landau-Zener and 
Kibble- Zurek), one can study the dynamics of large 
systems using time-evolving block decimation (TEBD) 
methods that use a matrix product state representation 
of the wavefunction. However, this approach is limited 
to ID, so we have not explored it further. 

Generalizations 

There is a great deal of flexibility in the implementa- 
tion of the corral idea. The corral can be formed either 
from a Bose Mott insulator, as described in the text, or 
from a Fermi insulator with one fermion localized at ev- 
ery site. Similarly, the corralled species can be either a 
Bose MI or a Fermi BI (as described). 
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FIG. 3: Time evolution of mean energy -Eavg, and energy uncertainty -Ewid, staggered structure factor Sq, and fidelity / as a 
function of time-dependent coupling J', for various ramp times tramp J = 0.25, 1, 3, 10, 30 (colors), on lattices of various sizes. 
The lattice sizes, from top to bottom, are 4 x 2, 4 x 4, 4 x 6, and 6x4 (the last of which was presented in the main text). 
Dashed lines on Sq plots indicate N^i^_^g/12, the value expected for a classical Neel antiferromagnet with staggered moment 
averaged over a sphere. On small lattices the true Sq can exceed this classical value due to the peculiar nature of quantum 
mechanical spin. 




FIG. 4: An alternative site-doubling procedure, in which a 
sublattice is gradually ramped up to introduce sites in the 
centers of the original unit cell. 



